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ABSTRACT 


This study investigates the mechanical characteristics of aluminum alloy 
AA3003-H14 when it is subjected to dynamic axial loads at low strain rates of less than 
1 s'. The tensile experiments show that both the ultimate and yield stresses, as well as the 
absorbed strain energy, increase with strain rate. Moreover, the strains at yield, ultimate 
load and fracture show positive sensitivity to strain rate. On the other hand, when the 
material is subjected to dynamic loading rates of more than 10°! s", the elastic modulus 
diminishes as the strain rate increases, while both normalized yield and ultimate stresses 
increase noticeably. It is also appreciated that the amount of plastic strain energy 
accumulated in a ductile material in its loading time history is useful as a failure criterion 
for the prediction of failure. Finally, a failure criterion is proposed for loading with 


varying strain rates. 
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I. INTRODUCTION 


The behavior of materials under high rates of strain may be vastly differently 
from that under static loading conditions. As the strain rate increases, the failure mode 
may change and the material strength may vary. This behavior, which differs for different 
materials, remains to be thoroughly investigated and appreciated. While the static 
strengths of composite materials, for example, have been considerably explored and 


documented, fewer studies have investigated their performance under dynamic loadings. 


Research on the effect of strain rate on the strength of materials is important to the 
design of their applications. An unmanned aerial vehicle, for example, may experience 
significantly higher loads only during the launch and recovery phases of their flights. 
Since these are not sustained static loads, the design of such vehicles to meet strength 
requirements derived from treating these loads under static conditions may be excessively 


conservative. Higher weights result and performance specifications are thus lowered. 


A. LITERATURE REVIEW 


1, Dynamic Loading on Metals 


Hadianfard et al. [1] investigated the effect of the rate of strain on the mechanical 
properties and failure mechanisms of the aluminum alloys, AA5754 and AAS5182. 
Conducting quasi-static tensile tests at strain rates of less than 10° s’, their results 
indicated negative strain rate sensitivity in these alloys. Both ultimate and yield strengths 
were reduced as the strain rate was increased. Moreover, the flow stresses and the strains 
to failure tended to decrease with higher quasi-static strain rates. The serrated yielding 
phenomenon associated with the nucleation and propagation of deformation band along 
the tensile specimen was observed. On the other hand, at dynamic rates of strain, the 
studied alloys exhibited mild positive sensitivity. The elongation to failure increased with 


the strain rate. 


Accordingly, different failure mechanisms were also observed in the quasi-static 
and dynamic tensile tests. Under quasi-static conditions, strain localization and shear 
band formation were found to be necessary pre-requisites for damage and final failure to 
occur. The thickness of the shear bands, as well as the average size of the damaged 
particles, was observed to increase with increasing strain rates. The contrary was true 
under dynamic conditions. The void nucleation, growth and coalescence process at the 
second phase particles was the dominant failure mechanism under dynamic conditions. 


Shear banding was less pronounced. 


Hadianfard et al. attributed the negative strain rate sensitivity in the quasi-static 
tensile tests to the dynamic strain ageing phenomenon which was manifested in the 
serrations. Being diffusion controlled, dynamic strain ageing was influenced by 
temperature and strain rate. Higher flow stresses at very low strain rates were the 
consequence of solute atoms interacting with obstacles preventing dislocation 
movements. With increasing strain rate, this dampening effect was thus reduced and 


lower flow stresses resulted. 


Similar results were also obtained by Mukai et al. [2] in fine-grained IN905XL 
aluminum alloys when they examined the dependence of their mechanical properties on 
the rate of deformation. When dealing with up to a strain rate of 10 s’, all samples 
showed negative strain rate sensitivity of flow stress. Positive strain rate sensitivity of 
strength was, however, observed above strain rates of 10° s!. In addition, while yield 
stress and total elongation were weak functions of the strain rate below a strain rate of 10° 
' |, both increased significantly with higher strain rates when the strain rate exceeded 
10° s7, 

In their tensile experiments on the aluminum alloy 6061, Srivatsan et al. [3] 
varied the strain rate from 10“ s" to 107 s”'. Likewise, they observed negative strain rate 
sensitivity, with yield and ultimate strengths diminishing as strain rate increased. 
Concomitantly, both the elongation to failure and the reduction in area increased. These 
observations were the result of shorter test times at high strain rates. Hence, interactions 
between dislocations and the primary hardening precipitates (Mg2Si) and interactions 


between dislocations and the coarse constituent particles were limited. The material 
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strength was thus lower and this, together with fewer micro-cracks due to reduced 
interactions between dislocations and inter-metallic particles, resulted in enhancements in 


the elongation to failure. 


Yu and Jones [4] conducted tensile tests on both aluminum alloy and mild steel 
specimens. They carried out quasi-static experiments at strain rates of less than 2 x 10° s 
' and dynamic tests up to a maximum strain rate of 140 s'. Their observations showed 
that the flow stress of aluminum alloy was insensitive to the rate of deformation. The 
latter affected only the rupture conditions, augmenting both the engineering and true 


rupture strains with increasing strain rate. 


Results on the mild steel specimens indicated that strain rate had a weaker effect 
on mild steel at large plastic strains than at small plastic strains. In addition, although the 
true rupture strain was nearly unaffected by strain rate, the true rupture stress increased 
about 20% when the strain rate was increased from 10° s"' to 10° s’. Good agreement 
was found when the Cowper-Symonds constitutive equation which is reproduced as 


Equation (1), was fitted to the experimental data. 
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pug 
o=o0,(e)[I+—]’ (1) 
D 
60(&) refers to the static stress-strain relation while D and p are material constants. 


The dynamic uniaxial stress-strain relationship was also studied by Albertini and 
Montagnani [5] in austenitic stainless steels. Ranging the strain rate between 107 s and 
5 x 10° s', they found that higher strain rates led to increasing flow stresses but 


decreasing uniform and fracture elongations. 
Zz Dynamic Loading on Composites 


Armenakas and Sciammarella [6] reported experimental findings on the 


mechanical properties of glass fiber reinforced epoxy plates subjected to high rates of 


strain in the direction of the fibers. Analyzing their results, they established that the 


dynamic elastic modulus varies linearly with the logarithm of the strain rate. 


Dynamic loading tensile tests were also performed by Lifshitz [7] on balanced 
angle ply glass fiber reinforced epoxy composites. Failure stresses were observed to be 
considerably higher than the corresponding static values whereas failure strains and 


moduli were insensitive to the rate of loading. 


Daniel et al. [8] tested and characterized unidirectional graphite/epoxy laminates 
in the form of thin ring specimens at high strain rates. The dynamic modulus exceeded 
the static modulus by approximately 20 percent while the dynamic strength and ultimate 


strain were close to the static properties. 


Applying a novel technique to the tensile impact testing of both uniaxially 
reinforced carbon fiber/epoxy and woven-roving reinforced glass fiber/epoxy composites, 
Harding and Welsh [9] showed that the modulus, the fracture strength and the failure 
mode of the carbon fiber/epoxy composite were dependent of the strain rate. In contrast, 
impact rates of strain lead to significant increases in the failure strength, the failure strain, 
the absorbed failure energy and the modulus of the glass fiber/epoxy composite. Unlike 
the carbon fiber/epoxy composite in which damage was confined to areas near the 
fracture plane, damage in the glass fiber/epoxy composite extended over the gauge region 
with increasing strain rates. Extensive debonding also occurred between the fibers and 


the matrix. 


In their literature survey on dynamically loaded composites, Melin and Asp [10] 
noted the difficulty in comparing results from various investigations since the assessed 
strain-rate dependencies were prone to influence from variations in fiber volume fractions 
and materials. Their research investigated the effects of the strain rate on the transverse 
tension properties of a carbon fibers/epoxy composite. By varying the strain rate between 
100 and 800 s"', they found weak or no dependence of the transverse mechanical 


properties on the strain rate. 


Through tensile tests performed on glass/epoxy laminates at over a range of strain 


rates, Okoli and Smith [11] demonstrated that the Poisson’s ratio is insensitive to the 
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strain rate. This could be attributed to the presence of fibers in the composites. They also 
commented that the Eyring theory of viscosity described the strain rate effects on most 
unfilled polymers. The theory assumed that the deformation of a polymer involved the 
motion of a chain molecule over potential energy barriers and further suggested that the 


yield stress is a linear function of the logarithm of the strain rate. 


In another study on the effects of the strain rate, Okoli [12] carried out tensile, 
shear and three-point bend tests on a woven glass/epoxy laminate. A linear relationship 
between expended energy and the logarithm of the strain rate was indicated. This could 
be explained by the failure modes of the composite laminates. Since matrix yielding 
increased with strain rate, a greater role is played by the matrix in the fracture process 
and hence, more energy is expended. Results also suggested that both the shear and 


flexural energies to yield of the woven laminates varied linearly with the strain rate. 


As noted by Taniguchi et al. [13], most composite materials displayed strain rate 
effects on mechanical properties under loading in matrix-dominant directions. In 
addition, the strain rate dependence of the tensile strength was higher in specimens with 
higher fiber orientation. Their investigations showed the carbon fiber reinforced 
thermoplastic epoxy specimen fractured without a necking process under dynamic 
loading. This contrasted with quasi-static loading, under which the necking process 
resulted in a large fracture strain. Higher strain rates were observed to yield larger tensile 
strengths as well as elastic moduli, including transverse and shear. The Poisson’s ratio, on 


the other hand, decreased with the strain rate. 


Shokrieh and Omidi [14] examined the behavior of unidirectional glass fiber 
reinforced polymeric composites under uniaxial loading at quasi-static and intermediate 
strain rates of between 0.001 and 100s’. They reported that the tensile strength and the 
absorbed failure energy of the composites enhanced significantly with increasing strain 
rate. Moreover, the tensile modulus and the strain to failure both showed slight increases 
when strain rates are raised. Failure modes were also observed to change from quasi- 


static to high dynamic loading conditions. 





























Table 1. | Review of experimental results on the effect of strain rate on the 
mechanical properties of composites. 
Effect of increasing strain rate on 
: Ultimate 
Study by Composite Failure Other 
tensile Modulus ; ee 
strain characteristics 
strength 
Armenakas 
and 
; Glass/epoxy Decrease Increase Decrease - 
Sciammarella 
[6] 
ae Angle ply 
Lifshitz [7] Increase | Independent | Independent - 
glass/epoxy 
Daniel et al. Increase 
Carbon/epoxy | Independent Independent - 
[8] slightly 
Absorbed 
Harding and | Glass/epoxy Increase Increase Increase | failure energy 
Welsh [9] increases 
Carbon/epoxy | Independent | Independent - - 
Melin and Transverse Increase Increase 
; Independent - 
Asp [10] carbon/epoxy slightly slightly 
Poisson’s 
Okoli and ; 
Glass/epoxy Increase Increase - ratio 
Smith [11] : 
independent 
Absorbed 
Okoli [12] Glass/epoxy - - - failure energy 

















increases 




















Effect of increasing strain rate on 
: Ultimate 
Study Composite Failure Other 
tensile Modulus ; Sut 
strain characteristics 
strength 
Carbon/ Poisson’s 
Taniguchi et Increase 
thermoplastic ‘ Increase Decrease ratio 
al. [13] slightly 
epoxy decreases 
Absorbed 
Shokrieh and Increase Increase ; 
-_ Glass/epoxy Increase failure energy 
Omidi [14] slightly slightly : 
increases 




















It is noted that several of the above studies suggested linear relationships between 
the logarithm of strain rate and various material characteristics including dynamic elastic 


modulus, yield stress and expended energy. 


E =aloget+b (2) 
o, =cloget+d (3) 
w=eloget f (4) 


Equations (2) and (4) demonstrate the relationships in a general form. A more 


useful form of the equations is presented in Equations (5) to (7). 








i= Vies £. (5) 
Ei Zs 

0 
ae —1=/log Hs (6) 
y,0 Eo 
” _1=mlog = (7) 
Wy é;y 


In these equations, the dynamic elastic modulus, the yield stress and the expended 
energy are normalized against corresponding reference values evaluated or determined at 


a reference strain rate. The reference values and the reference strain rate are denoted by 
the subscript ‘0’. 
The Cowper-Symonds constitutive equation of Equation (1) may be rearranged as 


follows: 





Oo 1 € 
lo 1/=—lo 8 
{2 P BD (8) 


It may thus be seen that the main difference between Equations (6) and (8) is in 
the left-hand sides of the equations. In Equation (8), the logarithm of the normalized 
stress is taken but not in Equation (6). 

Taken together, Equations (2) to (4) imply a constant yield strain as the strain rate 


varies, as shown below. 


2 
1 [coger a 


aloge+b 





2 
i oe [tog é) +2cd log e+d* 


2 


aloge+b 
DP 2abcd —b*c? 
co - 2ed be? a 
=~} — log e+ — -— + ——™—§— 
i ss alogét+b 


Thus, by comparing coefficients with Equation (4), Equations (9) to (10) may be 





obtained. 

C 

2cd_ be? 

ee a 
op 28,2 

Pi oe b°c _¢ (11) 


In particular, Equation (11) may be further manipulated to give Equation (12). 


ad’ —2abcd +b’c? =0 


(ad —be) =0 
ad = bc 
pe 
Thus, 
a, 
o =— 
é,=— a 
“  E 
_ cloge+d 
alog é+b 
ra 
=—+ a 
aloget+b 
ed 
a 
= constant 


(12) 


Thus, in materials where yield strains are not constant with varying strain rates, 


this may be attributed to at least one of the following reasons over the relevant range of 


strain rates: 


(1) The elastic modulus is not a linear function of the strain rate in this 


material; 


(11) The yield stress is not a linear function of the strain rate in this material; 


(iii) | The strain energy is not a linear function of the strain rate in this material. 


B. OBJECTIVES 


An objective of the current study is to examine the mechanical characteristics of a 
material when it is subjected to dynamic uniaxial loads. It has been noted that few studies 
have reported the variation of the elastic modulus with strain rate in metals. Thus, one 


purpose of the current study is to ascertain this relationship through tensile experiments. 


With the appreciation of the dynamic loading characteristics of a material, the 
study seeks to establish the latter’s failure envelopes in term of strain rates and 
normalized failure strengths. In addition, this study also investigates the failure criterion 
when a material is dynamically loaded to low strain rates of less than 1 s'. To the 
author’s best knowledge, there has been no failure criterion proposed for a varying strain- 


rate loading condition. 


Chapter II of this report describes the experimental setup and procedures 
employed in this study. Results and discussion are contained in Chapter III while Chapter 


IV closes with conclusion and recommendations. 
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Hl. EXPERIMENTAL SETUP AND PROCEDURES 


The material of the specimens tested is aluminum alloy AA3003-H14. The 
dimensions of the specimens are given in Figure 1. The thickness of the specimens is 


3.175 mm. 


125 





ALL FILLET R 3 MIN 


Figure 1. | Dimensions of test specimens. All dimensions in millimeters. 


The dimensions of the test sections of the specimens are first measured using a 
pair of vernier calipers. Tensile experiments are thereafter carried out on an Instron 4507 
Universal Materials Testing Machine at room temperature with a 20 KN load cell. The 
results of the tests are retrieved using a Series IX Instron software. Two important 


limitations of the current experiments are a maximum crosshead speed of 500 mm/min on 
’ hueos 50 
the testing machine and a maximum data acquisition rate of = Hz. These restrict the 


maximum strain rate that may be tested in the experiments. 


In the first set of experiments, strain rate is varied between 10° s‘' and 5x 10's". 
Strain rate may be calculated from the crosshead speed of the testing machine using 


Equation (13). 


e=— (13) 


V and lo denote the crosshead speed of the testing machine and the gauge length 
of the specimen respectively. Thus, with a gauge length of 15 mm, crosshead speeds of 
0.9 mm/min, 9 mm/min, 45 mm/min, 90 mm/min, 270 mm/min and 450 mm/min are 
tested. Crosshead speeds and hence, strain rates, are maintained constant throughout each 


test. 


All experiments are repeated and consistency is verified between the two sets of 


data. Parameters are thereafter computed as the average of the two. 


Figure 2 shows a typical stress-strain curve obtained in these experiments. 








Sa Necking 


one 
| i 
™~ 


a 


0 fy = ut © fracture 


Figure 2.__ A typical stress-strain curve. 


A straight line passing through the origin is first fitted to the linear portion of the 
curve. The elastic modulus is thus obtained as the slope of the straight line. In 


determining the yield strength and strain, the offset method is used. A second straight line 
14 


parallel to the first is offset 0.2% to the right of the latter. The 0.2% yield point is hence 


located where the second line interests the stress-strain curve. Numerical integration is 


carried out to compute the area under the curve so as to evaluate the strain energy. 


In addition to the first set of experiments which are performed with constant strain 


rates in this study, a second set of experiments are carried out, in which the strain rate is 


altered in each test. More specifically, in these experiments, the test is started with an 


initial crosshead speed. This is kept constant before it undergoes a step change in 


crosshead speed to a second speed. The latter is again maintained constant until the 


specimen fails. Table 2 tabulates the corresponding strain rates and the strains at which 


the changes in strain rate occur in these experiments. 




















Table 2. Parameters in varying strain rate experiments. 
No. First strain rate (1/s) Second strain rate (1/s) Transition strain 
1 0.01 0.30 0.00556 
2 0.01 0.30 0.00592 
3 0.01 0.30 0.01059 
4 0.01 0.30 0.01023 
5 0.30 0.01 0.01930 
6 0.30 0.01 0.02132 














These varying strain rate experiments object to examine the criterion for material 


ductile failure under varying strain rate loading conditions. 
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HiIl. RESULTS AND DISCUSSION 


The variations of the ultimate and yield stresses with strain rate are shown in 
Figure 3. Both the ultimate and yield stresses increase at higher strain rates. In addition, 
the rates of increase are noticeably higher at lower strain rates. 

184 - 
182+ 


180 + 


178 |- + 
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Figure 3. Variation of the ultimate and yield strengths with strain rate. 
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Experimental data is fitted both to Equations (6) and (14) as well as to Equations 
(15) and (16). While agreement is good with either set of equations, better agreement is 
found with Equations (15) and (16). In addition, a more apparent difference between the 


two sets of equations will be shown later. 





Cut = plog a +1 (14) 
O1t,0 e 
% 4 
é ; 
Oo, = oa)" (15) 
- 4 
== é Putt 16 
Ou = Sy ol 1 +—] (16) 


The coefficients, 1 and n, are determined as 0.0112 and 0.0137 respectively. Dy, 
Py, Dutt and pur are evaluated as 885 s”, 2.24, 2260 s” and 2.68 respectively. 


The strains at yield, ultimate load and fracture show increasing trends with strain 
rate, as shown in Figure 4. It may be observed that the increases in fracture strain are 
more pronounced at strain rates of below approximately 107 s'. On the other hand, the 
yield strain remains relatively constant at strain rates of below approximately 10° s" and 


increases thereafter with strain rate. 
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Figure 4. Variation with strain rate of the strains at fracture, ultimate load and yield. 


The elastic modulus is plotted against strain rate on a logarithmic scale in Figure 5 
for a better illustration of its variation with strain rate. As may be appreciated from the 
figure, the elastic modulus remains relatively constant at strain rates of below 10" s". 
Above this strain rate, the elastic modulus follows an approximately linear negative trend 


with the logarithm of the strain rate. 


Figures 6 and 7 plot normalized yield and ultimate stresses respectively against 
the inverse of the strain rate. Results from various studies are compared with results from 
the current study. The stresses are normalized by a corresponding reference value in the 
particular study. This is taken as the yield or ultimate stress obtained in the test of the 


lowest strain rate. 
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Figure 5. Variation of the elastic modulus with strain rate. 


At low deformation rates, the data shows relative insensitivity of both the 
normalized yield and ultimate stresses to the strain rate. On the other hand, as the strain 
rate is increased beyond 10" s"', both the normalized yield and ultimate stresses increase 


noticeably. 


The inverse of the strain rate, which is plotted on the x-axis in the figures, gives 
an indication of the duration of time that the material is allowed to deform. Plastic 
deformation by loading beyond the yield strength or failure at ultimate strength, 
respectively, occurs in the material when the stress level lies in the shaded zones above 
and to the right of the plotted points. The plots show that when deformation is rapid and 
the material is not allowed a sufficient time to deform, plastic deformation and failure in 
the materials are delayed to higher stresses. This effect, however, may be observed to 


diminish when the strain rates are less than 107! s7. 
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Figure 6. Variation of normalized yield strength with the inverse of the strain rate. 


The dotted and dashed lines in the graphs represent extrapolated curves fitted to 
Equations (6), (14), (15) and (16). Generally, Equations (15) and (16) are better 
representations of the experimental data obtained than Equations (6) and (14). The 
dashed lines in Figures 6 and 7 thus define envelopes for plastic deformation and failure 


in terms of the strain rate and the normalized yield and ultimate strengths. 
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Figure 7. Variation of normalized ultimate strength with the inverse of the strain 
rate. 


The total strain energy density absorbed by the material to failure is plotted as a 
function of strain rate in Figure 8. It is shown to increase with increasing strain rates. 
However, the rate of increase diminishes at higher strain rates. Equation (7) may be fitted 
to the experimental data with good agreement, as may be observed in Figure 8. The 


coefficient, m, is obtained as 0.106. 


22 





L 


3.6 L L L L L L 
0 0.05 «80.1 0.15 O02 025 O83 O35 04 045 0.5 
Strain rate (1/s) 





(a) Strain rate in linear scale. 


4.8+ 


4.67 





3:6 3 3 1 o 
10 10 10 10 


Strain rate (1/s) 





(b) Strain rate in logarithmic scale. 


Figure 8. __ Variation of the total strain energy density with strain rate. 
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Total strain energy density comprises elastic and plastic strain energy densities. 
Figure 9 shows the effect of strain rate on the elastic strain energy density. As suggested 
in Figure 9, the elastic strain energy density may be expressed as a positive linear 


function of the strain rate. This is given in Equation (17). 


ere (17) 


The values of q and r are 0.433 and 0.0772 respectively. 
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Figure 9. ‘Variation of the elastic strain energy with strain rate. Inset shows the 
definition of elastic strain energy. 


The variation of the plastic strain energy density with strain rate is illustrated in 
Figure 10. Strain rate may be observed to have a similar effect on the plastic strain energy 


density as the total strain energy. The plastic strain energy density increases as the strain 
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rate increases. However, the rate of increase diminishes at higher strain rates. When 
Equation (18) is fitted to the experimental data, reasonable agreement is found with the 


coefficient, s, evaluated as 0.0966. 





—l=slog Bae (18) 


It is noted, however, that Equation (18) does not reconcile with Equations (7) and 


(17). This may be explained in the following. 


w, =w-w, (19) 


Equation (18) may be used to evaluate the left hand side of Equation (19) as 


below. 


LHS of Equation (19) =w, slog zt +w 
€ 


p.0 


=s'log = +w,, where s'=w, \s (20) 
€ 


On the other hand, when the right hand side of Equation (19) is evaluated with 


Equations (7) and (17), another expression for the plastic strain energy may be obtained. 
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Figure 10. Variation of the plastic strain energy with strain rate. Insets show the 
definition of plastic strain energy. 
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RHS of Equation (19) =w,mlog +W, gé r 
é 





! é : : ! 
=m'log| — |+w,—-qgé+qé,—w,, where m'=w ym 
& 
= fe 
=m'log] — |+w,)—-q(é- &)) (21) 


The difference between the above two expressions of the plastic strain energy 


density, that is, the difference between the left and right hand sides of Equation (19), lies 


in the last term of the right hand side, ~q(€-&). Due to setup limitations on the 


achievable range of strain rates, however, insufficient data is available in this study to 


resolve with confidence the validity of either expression. 


As a failure criterion in ductile materials, the use of a critical value of the plastic 
strain energy density is appropriate. Consider, for example, a metal repeatedly loaded to 
different stresses below the ultimate strength and unloaded thereafter. While elastic strain 
energy is recovered upon each unloading, the plastic strain energy represents a reduction 
in the strength of the material. Material failure occurs when its plastic strain energy 
density exceeds a critical value. In more general terms, Li [15] concluded that material 
ductile failure is influenced by two dissipative mechanisms, namely, the plastic 


deformation dissipation and the damage dissipation. 


In order to demonstrate the application of critical values of the plastic strain 
energy density as a ductile failure criterion, consider a general strain time history in 
Figure 11. In this strain time history, loading and unloading occurs at varying strain rates. 


However, for simplification, strain rate is discretized such that the strain rate of each 
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loading and unloading phase is approximated by the slope of a straight line drawn from 


the point of the lowest strain to the point of the highest strain. This is illustrated in Figure 


Ld. 


7 & 878 _ & ~ 5 Sk, 2-2-6 


oe 
In this example, consider a sequential strain rate loading such as ¢,=—, 
t, 


. €,—€ 
, 6 =-4 SG=4 t=+ and ¢,=4 2 
t-t; t-t t,t, t,-t, t.-t, t-te t, —t, 





where E > Ey > é, > E, . The time variation of these strain rates are plotted in Figure 12. 


Figure 13 shows the corresponding stress-strain curves due to the respective strain 


rates of E; , Ey : &; and E, . These may be determined using Equation (1) from the static 


stress-strain relation once the material constants are known. The variations of the stress 


and strain in the material with time resulting from the strain time history shown in Figure 


11 is illustrated in Figure 13 in blue. 





Figure 11. General strain-time history. 
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Figure 12. Time variation of strain rate due to strain time history. 
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Figure 13. Stress-strain curves due to strain time history. 
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In each loading cycle in which the material is loaded and subsequently unloaded, 
its plastic strain energy gained in that loading cycle is determined as the respective 
shaded area. In the computation of the plastic strain energy density, it is important to note 
that yield strain and elastic modulus are functions of strain rate. The plastic strain energy 
accumulates over the loading cycles and material failure is predicted when this energy 
density exceeds a critical value. In particular, a failure criterion for varying strain rate 


dynamic loading conditions is proposed in Equation (22). 





(w,), (w,), (w,), (w, ) 
= + +..+ n—] 22 
Gh Utah Oral. Wel 22) 


(w,) denotes the amount of plastic strain energy density stored in the material 


during the ith loading cycle while (w refers to the critical plastic strain energy 


p.critical ) : 
density at the same strain rate of the loading cycle. With the discretization of the strain 
rate, each loading phase is deemed as occurring with a constant strain rate. Hence, the 
critical plastic strain energy density is a function of the strain rate of the loading phase 


and may be given either by Equation (20) or Equation (21). 


The proposed failure criterion expressed in Equation (22) is corroborated by 
results obtained in the varying strain rate experiments. Figure 14 shows a typical stress- 
strain curve obtained in such an experiment. In this instance, the strain rate is increased 
from 10° s' in the first stage of the experiment to 3 x 10°! s” in the second stage at a 
transition strain of 0.01023. The stress-strain curves due to the tensile tests of constant 
strain rates of 10° s' and 3 x 10° s are shown in Figure 14 for comparison. As 
expected, the stress-strain curve follows that of the 107 s” strain rate initially. Upon the 
step increase in strain rate to 3 x 10's", the stress-strain curve begins to deviate from 
that of the 107 s"' strain rate to lie in between the two stress-strain curves of constant 


strain rate. 
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Figure 14. _A typical stress-strain curve obtained in varying strain rate experiment. 


In efforts to assess the validity of the proposed failure criterion, the plastic strain 


energy densities absorbed by the material during the two stages of each varying strain 


rate test, (w,), and (w, 


D 


iz are computed. These are determined as the shaded areas as 


illustrated in Figure 14 using numerical integration. Each stored plastic strain energy 
density is divided by its respective critical plastic strain energy density, and the two ratios 
are summed. Results are tabulated in Table 3, which shows that the proposed failure 
criterion in Equation (22) is validated by these experimental results within an error of 


approximately 10%. 


In another approach, the proposed failure criterion is applied to each of the 
varying strain rate experiments to predict fracture strain. These values are estimated 
based on the transition strain as well as the constant strain rate stress-strain curves of the 


first and second strain rates. It is not presumed that stress and strain measurements are 


a1 


made in the varying strain rate experiments. Table 4 compares the computed results with 


the empirical data, and this comparison indicates that the former approximates the latter 


within errors of approximately 4%. 






















































































Table 3. Results from varying strain rate experiments. 
Second 
First strain ; Transition (w, ), (w, (w, ) 
No. strain rate cma), 
rate (1/s) (1/s) strain (MJ/m?) (MJ /m) Wo critica! J, 
1 0.01 0.30 0.00556 0.5482 3.6273 0.9461 
2 0.01 0.30 0.00592 0.6540 J20238 0.9501 
3 0.01 0.30 0.01059 1.4644 2.8404 1.0057 
4 0.01 0.30 0.01023 1.4130 2.7365 0.9695 
5 0.30 0.01 0.01930 2.7684 1.1901 0.9196 
6 0.30 0.01 0.02132 3.3526 0.6405 0.9086 
Table 4. Actual and predicted fracture strains in varying strain rate experiments. 
Second Actual Predicted 
First strain Transition Percentage 
No. strain rate ; fracture fracture ; 
rate (1/s) strain difference (%) 
(1/s) strain strain 

1 0.01 0.30 0.00556 0.0270 0.0279 3.1681 
2 0.01 0.30 0.00592 0.0268 0.0278 3.7789 
3 0.01 0.30 0.01059 0.0279 0.0272 -2.5958 
4 0.01 0.30 0.01023 0.0268 0.0272 1.6833 
5 0.30 0.01 0.01930 0.0274 0.0268 -2.1726 
6 0.30 0.01 0.02132 0.0262 0.0269 2.6252 
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More specifically, the approximation methodology is described in the following. 
Firstly, the stress-strain curve for the constant strain rate of the first stage is integrated 


from zero to the transition strain so as to obtain the plastic strain energy density gained 


(w, ), 
(Was ) 


the stress-strain curve for the constant strain rate of the second stage is integrated from 


during the first stage. is thereby checked to be less than unity. Subsequently, 


the transition strain to an upper limit that is continually increased until the failure 


criterion is satisfied. The fracture strain is thus estimated by the upper limit. 


In the case of the last two experiments in Table 4, when the stress-strain curve for 


the constant strain rate of the second stage is integrated from the transition strain to its 


P P 
1 2 
( W > critical ), ( Wp critical ), 


to predict the fracture strain, a fourth order polynomial is fitted to the plastic deformation 





fracture strain, remains less than unity. Thus, in order to be able 


regime of the strain-strain curve for the constant strain rate of 10 s' and extrapolated. 
This allows integration to be carried beyond the fracture strain due to the constant strain 


rate of 107! 57. 


Additional experiments tend to suggest that the strain rate of unloading influences 
the variation of stress with strain during an unloading. When a material is unloaded at a 
particular strain rate, its stress varies as a linear function of strain. The slope of this 
function is perceived to be the same as the elastic modulus of that strain rate. Thus, if a 
material is loaded at a different strain rate than when it is unloaded, the slope during 


elastic deformation would differ from the slope during unloading. 


The range of the strain rates tested in this study spans two different regimes. The 
quasi-static and dynamic regimes are classified, respectively, with strain rates of less than 
and more than 10°! s”. This is similar to observations by Hadianfard et al. [1] who also 
differentiated between tensile tests at strain rates of less than and more than 10°! s"'. In 


the current study, differences between the two regimes are evident in Figures 4 to 7. 
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There are two limitations in the above result analysis. Firstly, due to a lack of data 
on true stress and strain, numerical integration of the area under the engineering stress- 
strain curve is carried out for the computation of strain energy density. Secondly, due to 
constraints on the maximum rate of data acquisition, the number of data points acquired 
in the tests of the highest strain rate is small. In particular, this reduces the accuracy of 
the calculated elastic modulus, yield stress and yield strain in the tests of strain rate 5 x 
10's". 

For completeness, Figure 14 is presented as the corresponding graph to Figure 7 
in composite materials. Similar trends may be observed. When deformation is rapid and 
the material is not allowed a sufficient time to deform, material damage is delayed to 
higher loads. This effect, however, may be observed to diminish at strain rates of more 


than 107s 
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Figure 15. Normalized ultimate strength as a function of the inverse of the strain rate. 


34 


IV. CONCLUSIONS AND RECOMMENDATIONS 


In conclusion, the mechanical characteristics of aluminum alloy AA3003-H14 
were examined under a dynamic uniaxial loading. Both ultimate and yield strengths, as 
well as the absorbed strain energy density, show positive sensitivity to strain rate. 
Moreover, the strains at yield, ultimate load and fracture increase with strain rate. Under 
quasi-static strain rates of less than 10°! s’, the low loading rate does not have a 
significant effect on the elastic modulus, the normalized yield stress or ultimate stress. 
On the other hand, when the material is subjected to dynamic loading rates of more than 
10's", the elastic modulus diminishes as the strain rate increases, while both normalized 
yield and ultimate stresses increase noticeably. It was also noted that the amount of 
plastic strain energy density accumulated in a ductile material in its loading time history 
was useful as a failure criterion for the prediction of failure. As a result, a failure 
criterion was proposed for loading with varying strain rates using the concept of 
accumulated plastic strain energy density. The proposed criterion was validated from 
experiments consisting of two different strain rates applied one after another 


subsequently. 
The current work may be extended in the following ways: 


(1) The range of strain rates used in the tensile experiments may be increased. 
This would enable the verification of the trends observed in the current study beyond the 
present maximum strain rate. In particular, this would allow the resolution of an 


empirical expression for the plastic strain energy. 


(ii) Different metals may be investigated and tested. By studying the results 
due to different metals, the variations of the normalized yield and ultimate stresses with 
the inverse of the strain rate may be compared. Their similarity or differences may thus 


be established. 


(iii) | Composite materials may be investigated and tested. Where composites 
fail in a brittle manner, the effects of strain rate on their mechanical properties may differ. 
Moreover, a different failure criterion may apply. 
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(iv) Much remains unclear about the behavior of materials under dynamic 
loading and unloading conditions of varying strain rates. More experiments designed to 
investigate the variation of stress with strain under these conditions would enhance the 
understanding of the influence of varying strain rate on the evolution of plastic strain 


energy. 
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